function [C1] = Ddavid_estimate_sampling_constant_1(VLabel, VResult)

M = size(VLabel, 2);
C1 = zeros(1, size(VLabel, 2));

for i = 1:M
    NP = size(VLabel(:, i) == 1, 1);
    LabelSum = sum(((VResult(VLabel(:, i) == 1, i) + 1) / 2), 1);
    C1(1, i) = LabelSum / NP;
end
